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Motivated by recent Hall-effect experiment in YbRli2Si2, we study ground state properties of a 
Kondo lattice model in a two-dimensional square lattice using variational Monte Carlo method. 
We show that there are two types of phase transition, antiferromagnetic transition and topological 
one (Fermi surface reconstruction). In a wide region of parameters, these two transitions occur 
simultaneously without the breakdown of Kondo screening, accompanied by a discontinuous change 
of the Hall coefficient. This result is consistent with the experiment and gives a novel theoretical 
picture for the quantum critical point in heavy fermion systems. 
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Quantum phase transition where a transition temper- 
ature goes to zero is a very interesting phenomena often 
observed in strongly correlated electron systems and has 
attracted much interests these days. Compared with the 
conventional finite-temperature phase transition, quan- 
tum nature of fluctuation explicitly appears. In heavy 
fermion systems, the quantum phase transition occurs 
when the antiferromagnetic (AF) long-range order dis- 
appears, as some parameters such as pressure, doping 
and applied magnetic field changes. Around this quan- 
tum critical point (QCP), many interesting phenomena 
such as non-Fermi-liquid behavior or unconventional su- 
perconductivity have been observed. Recent experimen- 
tal development makes it possible to study the property 
of the QCP in detail and suggests the necessity of the 
improvement of a conventional view for the QCP. 

Recently, it was found that the Hall coefficient of a 
heavy fermion material, YbRfi2Si2, shows a rapid change 
as a function of magnetic field When the tempera- 
ture is lowered, this crossover point approaches the QCP 
and simultaneously the crossover becomes much sharper. 
From this result, it was claimed that the crossover be- 
comes a discontinuous jump at T = and that a re- 
construction of the Fermi surface (FS) occurs at the 
QCP, i.e., from a large FS to a small FS Q. In the 
conventional view of heavy fermion systems, the conduc- 
tion electrons (c-electrons) and the localized spins (/- 
electrons) hybridize with each other through the Kondo 
coupling and form "heavy Fermi liquid" on both sides 
of the QCP. In this case, the volume of the FS is deter- 
mined by the total number of the c- and /-electrons and 
thus the system has a "large FS" . Once such a Fermi liq- 
uid is formed, the self-consistent-renormalization (SCR) 
theory by Moriya Q and renormalization-group studies 
by Hertz and Millis [4] can be applied. However, the 
FS reconstrunction observed experimentally can not be 
understood. Recently Si proposed an alternative theory 
for the QCP, called as "local quantum criticality" [H, Q. 



He claimed that the Kondo screening becomes irrelevant 
at the QCP and the c- and /-electrons are decoupled in 
the AF phase. As a result, /-electrons turn to be local- 
ized and do not contribute to the FS any more. If this is 
the case, the volume of the FS is determined only by the 
number of c-electrons and a so-called "small FS" is real- 
ized. Although this picture of local quantum criticality 
seems reasonable to explain the Hall-effect experiment, 
it is still sketchy and the details of the theory are contro- 
versial. 

In this paper, we study this problem microscopically 
by applying a variational Monte Carlo (VMC) method 
to the Kondo lattice model (KLM). Since the Hall-effect 
experiment suggests a zero-temperature phase transition 
at QCP, the T = variational approach is suitable to this 
problem. The VMC method has been successfully used 
in the strongly correlated electron systems. We construct 
several trial wave functions which have large- or small FS 
together with AF long-range order. Evaluating the vari- 
ational energies with VMC simulation, we determine the 
ground state phase diagram as a function of the Kondo 
coupling J and the c-electron number n c . We find a first 
order phase transition from a "large FS" to a "small FS" 
as a function of J, even if the Kondo screening always 
occurs. The present results give a new route of the FS 
reconstrunction which was not considered in local quan- 
tum criticality. 

We study the following KLM in a two-dimensional 
square lattice, 



H 



(i) 



where Si and Si represent the /- and c-electron 
spins, namely, Si = | fi* a ™' f™' and s * = 
\ ^aa' c L°W c icr'- J(> 0) denotes the AF exchange 
coupling between them. As for the kinetic energy term, 
we consider the nearest-neighbor hopping and thus the 
bare c-electron energy dispersion is given by 



e& = — 2t(cos k x + cos k y ). 



(2) 
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In principle, both the Kondo effect and the RKKY in- 
teraction are the consequences of the "intrasite" inter- 
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(a)PM (b)AFh (c)AFe (d) AF S 

FIG. 1: Fermi surfaces and band dispersions of the four different variational states considered in this paper, (a) paramagnetic 
metal (PM), (b) AF h , (c) AF e and (d) AF S . 



action between c- and /-electron spins. However, many 
mean-field type theoretical studies include an additional 
term of the "intersite" exchange interaction, Jh Si-Sj, 
for theoretical convenience (This model is called "Kondo 
Heisenberg model" ). We have to deal with this additional 
term with care, since it sometimes leads to overestimate 
of the RKKY interaction and even to unphysical results. 
Therefore, we study the model without the intersite ex- 
change term in this letter. 

The KLM has been studied in a mean-field approx- 
imation (MFA) 0, 0. The disadvantage of MFA is 
that the local constraint of one electron per / orbital, 
^Zcrficrfia = 1, is not satisfied locally but only in aver- 
age. Indeed, other methods which can correctly enforce 
this local constraint at = 1.0 have shown that the 
MFA gives wrong results [H, [To[. In contrast, the VMC 
method takes account of the local constraint exactly. In 
the following, we focus on the case with n c < 1.0 where 
the issue of large- and small FS can be discussed. For 
this case (n c < 1.0), the conventional Quantum Monte 
Carlo method [§] can not be used due to the severe sign 
problem. The advantage of using the VMC method is 
twofold: (1) it has no negative sign problem, and thus 
can be applied to the system with n c < 1.0 or with geo- 
metrical frustration, and (2) the local constraint of each 
site is strictly satisfied. This constraint is necessary for 
the study of the KLM, but hard to be enforced by other 
methods. 

For the trial wave functions, we use the following 
Jastrow-type wave functions, 



where 



^=i=nK(Wi)+<(W T ) 



(3) 



(4) 



is a projection operator which keeps the /-electron num- 
ber of each site exactly one. |3>) is obtained by diago- 
nalizing a one-body Hamiltonian with some variational 
parameters. This type of wave function has been studied 
by Shiba and Fazekas for the one-dimensional case [ill ]. 
However, AF long-range order is not stabilized in one- 
dimension, so that the AF QCP was not studied [l2Lll3l|. 
In contrast, we show that the AF state is stabilized in a 
finite region of J < J c and we can discuss the AF QCP. 

First, let us explain the trial wave functions used in this 
paper in some detail. For a paramagnetic metal (PM) 
state, we construct |<1>) from the one-body Hamiltonian 



ff = E( c L./L) 
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where fka is a Fourier transform of fi a . V and Ef are 
variational parameters that control the degree of c-f hy- 
bridization and the effective /-electron level, respectively. 
Since the ratios V/t and Ef/t appear in the trial wave 
function, we do not need to treat t in ek as a variational 
parameter. We set t = 1 as a unit of energy in the fol- 
lowing. 

In the same way, AF state can be obtained from 

H = ( C Afccr> C Bfccr> /Ifccr' fuka 
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where m and M denote the variational parameters of 
AF moment in c- and /-electrons, respectively. A and B 
represent the indices of each sublattice. The summation 
over k is in the folded AF Brillouin zone. If we set m = 
M = 0, Eq. © reduces to Eq. ©. 
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FIG. 2: Ground state phase diagram in the J — n c plane ob- 
tained in VMC. Solid lines ( J c ,fs and J c ,af+fs) represent the 
first-order phase transition and dotted line ( J c ,af) represents 
the second-order one. 



We optimize the variational parameters mentioned 
above and find the lowest energy state in a J — n c phase 
diagram. Possible candidates for the ground state are 
classified according to the shape of the FS and the band 
dispersion, as shown in Fig. [TJ (a) paramagnetic metal 
(PM), (b) AF metal with hole-like FS (AF h ), (c) AF 
metal with electron-like FS (AF e ) and (d) AF metal 
without c-f hybridization (AFs, S denotes "small"). In 
PM, AFh and AF e states, V is finite, i.e., the c- and f- 
electrons hybridize with each other through the Kondo 
screening. The state called AFs represents the small FS 
proposed by Si, which has no Kondo screening, V = 0. 
The state PM is the conventional Fermi liquid state with 
large FS, where quasiparticles become heavy due to the 
hybridization and projection operator P / = 1. When 

AF long-range order is introduced into PM, we obtain 
AFh- Therefore, PM and AFh are smoothly connected 
with each other and their FSs have the same topology 
(hole-like FS). On the other hand, AF e (Fig. ffl(c)) is a 
state which can be obtained by introducing V in the AFs 
with small FS. This state have not been obtained before 
in the MFA, but it is a natural extension of AFs- The 
topology of FS of AF e is different from that of AFh as 
shown in Fig. U[b) and (c). However, since the Brillouin 
zone is halved due to the presence of AF long-range or- 
der, the FS volumes are the same. Actually, by tuning 
V, Ef, m and M, the band dispersion changes contin- 
uously from AFh to AF e (Fig. [1]). In this sense, there 
is no difference of symmetry between AFh and AF e , but 
the Hall coefficient will be different between these states. 
We regard the state AFh as a state with a "large FS" 
connected to PM (Fig. UK a)), and AF e as a state with 
"small FS" connected to AF S (Fig. [2(d)). 

Figure [2] shows the obtained phase diagram in the 
J — n c plane. The numbers of sites are set to be 64, 
100, 144 and 196 for the calculation and the results show 
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FIG. 3: Energy difference of each state compared with the 
PM state for (a) n c = 0.9375 and (b) n c = 0.8125. The error 
bars are order of 10 -4 and not shown here. 



good convergence. For n c > 0.82, we have two differ- 
ent types of transition. When J decreases from a large 
value, the variational ground state changes from PM to 
AFh at J c ,af- This is a "conventional" second-order AF 
transition where the AF order parameter grows contin- 
uously from zero to a finite value. The first Brillouin 
zone is folded and the number of energy band doubles 
(2 — > 4) as shown in Fig. [2(a) to (b). The energy differ- 
ence compared with the PM state is shown in Fig. (3](a). 
When this transition occurs, the Hall coefficient must 
change continuously. Then, as J decreases further, the 
state changes from AFh to AF e at J c ,fs- We find that 
this transition is first-order even if AFh and AF e are con- 
nected in the variational parameter space. Actually, the 
variational energy has a double minimum, and an energy 
crossing occurs as shown in Fig. [31(a) at J c± 1.1. In other 
words, the band dispersion changes discontinuously from 
convex downward to convex upward as shown in Fig. [2(b) 
and (c). It is regarded as a kind of Lifshitz transition, 
namely, the topology of FS changes. In this case, the 
discontinuous change of the Hall coefficient occurs. It 
is a novel type of quantum phase transition originated 
from the competition between the Kondo effect and the 
RKKY interaction. In addition to the Hall coefficient, 
the expectation value of each term in Eq. (pQ) shows a 
discontinuous change at J c ,fs- 

When n c < 0.82, AFh is not stabilized and the state 
changes directly from PM to AF e at J c ,af+fs- The ob- 
tained variational energies are shown in Fig. [31(a). In this 
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TABLE I: Several pictures for the QCP in heavy Fermion material YbRli2Si2. 



Picture 


Transition 


Kondo screening 


Change of Hall coefficient 


Experiment 


conventional QCP 


PM -> AF h 


remains 


continuous 


inconsistent 


local quantum criticality (Si et al.) 


PM -> AF S 


disappears 


discontinuous 


consistent 


our study (n c < 0.82) 


PM -> AF e 


remains 


discontinuous 


consistent 



case, the AF transition and the topological phase transi- 
tion of FS occur simultaneously. The AF order parameter 
and the Hall coefficient change discontinuously reflecting 
the character of the first order phase transition. 

Interestingly, there is no region where AF$ is stable. 
Namely, the Kondo screening always occurs in our phase 
diagram. This result is quite different from those ob- 
tained in MFA @, Hj] , where the Kondo screening is sup- 
pressed across the QCP and the AF$ becomes stable in a 
wide region. We think that this discrepancy appears from 
the insufficient treatment of the local constraint in MFA. 
As discussed before for the case with n c = 1.0 0,03, the 
AF order is overestimated in MFA and thus the MFA is 
not appropriate for the problem considered here. 

Finally let us compare the Hall-effect experiment with 
the obtained results. We summarize the theoretical pic- 
tures for the QCP in Table U In a conventional pic- 
ture, the Hall coefficient varies continuously across the 
QCP [14]. In our notation, the state changes from PM 
to AFh- In this case, even if there is a nesting of FS, the 
change is continuous unless the FS is perfectly fiat [151 ]. 
Therefore, this conventional picture is inconsistent with 
the experiment. On the other hand, we have found that 
the AF transition and the discontinuous change of the 
Hall coefficient occur simultaneously for n c < 0.82. This 
is consistent with the experiment, although it is not con- 
clusive whether the Hall coefficient has a discontinuity at 
T = experimentally. We would like to emphasize that 
our transition is different from the local quantum criti- 
cality by Si in the sense that the Kondo screening always 



occurs. We have shown that the FS topology changes dis- 
continuously even in the presence of hybridization. Our 
present results are T = variational calculation, and 
can not be applied directly to finite temperature. Hall 
effect experiment was carried out at finite temperatures 
and only a crossover takes place as a function of applied 
magnetic field. The application to finite temperature be- 
yond MFA is a challenging problem in future. 

In summary, we have studied the ground state of the 
KLM in a two-dimensional square lattice with the VMC 
method. We proposed new types of wave functions (AFh 
and AF e ) which have different FS topologies. It is shown 
that there are two kinds of phase transition, the AF 
transition and the topological one (FS reconstruction) 
even if the Kondo screening always occurs. We find that 
these two transitions occur simultaneously at the QCP 
for n c < 0.82 and this can explain the result of the Hall- 
effect experiment in YbRli2Si2. It is a novel type of quan- 
tum phase transition induced by the competition between 
the Kondo effect and the RKKY interaction. 
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